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The single-wave model equations are transformed to an exact hydrodynamic closure by using a class of 
solutions to the Vlasov equation corresponding to the waterbag model. The warm fluid dynamic equations 
are then manipulated by means of the renormalization group method. As a result, amplitude equations for 
the slowly varying wave amplitudes are derived. Since the characteristic equation for waves has in general 
three roots, two cases are examined. If all three roots of the characteristic equation are real, the amplitude 
equations for the eigenmodes represent a system of three coupled nonlinear equations. In the case, where the 
dispersion equation possesses one real and two complex conjugate roots, the amplitude equations take the 
form of two coupled equations with complex coefficients. The analytical results are then compared to the 
exact system dynamics obtained by solving the hydrodynamic equations numerically. 
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I. INTRODUCTION 

The processes of pattern and coherent structures for- 
mation in plasmas and plasma-like media have attracted 
attention for many years. A number of approxima- 
tions and simplified models has been proposed, amongst 
which the single-wave model is one of the most efficient 
approaches to study the weakly nonlinear behavior in 
marginally stable plasmas. Starting from the general 
Vlasov-Maxwell equations for the phase space density 
distribution and the self-consistent electromagnetic field, 
one usually derives a coupled set of equations describing 
the evolution of the coarse-grained distribution function 
and a certain isolated marginally stable wave mode of the 
electrostatic potential. 

Using the method of matched asymptotic expansions, 
the single-wave model equations have been recently 
derived^ in the most general case independent of the equi- 
librium. A rigorous derivation of the nonlinear evolution 
of a unstable electrostatic wave for a multispecies Vlasov 
plasma has been given by Crawford and Jayaraman^ and 
this general physical picture has been shown to corre- 
spond to the single-wave model. Amongst earlier pio- 
neering work utilizing the single-wave model, the stud- 
ies dedicated to the beam-plasma instability 3 ^- and the 
bump-on-tail instability 6-8 should be mentioned. The in- 
fluence of the finiteness of the number of particles coupled 
to a monochromatic wave in a collisionless plasma has 
been investigated by M-C. Firpo et al^ilil. Based on the 
utilization of the single-wave picture some analogies be- 
tween two classical models of the single-pass free electron 
laser dynamics and of the beam- wave plasma instability 
have been discussed by A. Antoniazzi et alii. This im- 
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plies that the results presented in the present paper will 
be useful to describe processes of formation of patterns 
and coherent structures in single-pass free electron lasers. 

The purpose of the present paper is not to argue on 
the validity or the justification of the single- wave picture, 
rather than to use it as a starting point for deriving am- 
plitude equations for the slowly varying wave envelopes. 
It is organized as follows. In the next section, we cast 
the self-consistent single-wave model equations derived 
by del-Castillo-Negrctc 1 in an equivalent form by using 
a class of exact solutions to the Vlasov equation corre- 
sponding to the waterbag mode l 12 ' 13 . Further, the exact 
closure of hydrodynamic equations are manipulated fol- 
lowing the renormalization group (RG) approachi 3 -. As 
a result, amplitude equations for the slowly varying wave 
envelopes (SVWE) are derived. Depending on the char- 
acter of the solutions of the dispersion equation, two cases 
can be distinguished. In the first case, where all three 
roots of the dispersion equation are real, the amplitude 
equations represent a system of three coupled nonlinear 
equations, describing the process of three-wave interac- 
tion and mixing. If the dispersion equation possesses 
one real and two complex conjugate roots, the amplitude 
equations take the form of two globally coupled nonlin- 
ear equations with complex coefficients. In Section V, the 
predictions of the SVWE equations are benchmarked to 
the exact numerical solution of the hydrodynamic equa- 
tions. Finally, in the last section we draw some conclu- 
sions and outlook. 



II. THEORETICAL MODEL AND BASIC EQUATIONS 

The basic equations derived by del-Castillo-Negretei, 
which will be the starting point of our subsequent anal- 
ysis can be written as 



dtf + VdJ + d x ipdvf = 0, 



if = a(t)e ix + a" ! (t)e~ ix , 



(1) 



(2) 
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a^+ila = i(fe- i % 
where the operator averaging is denoted by 

OO 2-7T 

1 

2^ 



(...) 



dV / dr . 





(3) 



(4) 



The independent time t and spatial x variables, as 
well as the dependent ones f(x, V;t) and <fi(x;t) 
entering the above equations have been properly 
nondimensionalizedi. For example, the dimensionless 
time t scales as u e T ', where uj e is the electron plasma 
frequency and T is the real time. Moreover, the velocity 
V in the equations above scales as V/Vd, where Vd is a 
velocity characteristic for the particular problem and can 
be chosen arbitrarily. 

For present purposes, we restrict the analysis to a spe- 
cial class of exact solutions to Eq. (JTJ corresponding to 
the waterbag distributio n 12 ' 13 

f{x, V; t) = C[U{V - v.(x; t)) - H(V - v+(x; t))}, (5) 

where T-L denotes the well-known Heaviside function, C is 
a normalization constant and < x < 2?t is a normalized 
spatial variable. It simply means that the phase space 
density /(x, V;t) remains constant within a region con- 
fined by the boundary curves v±(x; t), which are assumed 
to be single- valued. The latter distort nonlinear ly during 
the evolution of the system as specified by Eqs. (HJ - ©. 

It is convenient to introduce the macroscopic fluid vari- 
ables 



dVf{x,V;t) = C(v+-V-), 



(6) 



gv= I dVVf(x,V;t) = -(v 2 + -v 2 _), (7) 



where g(x;t) and v(x;t) are the density and the current 
velocity, respectively. It can be shown that the higher 
moments defining the particle pressure V{x;t) and the 
heat flow Q(x; t) can be expressed as 

OO 

V= J dV(V-vff(x,V;t) = ^(v + -v^f, (8) 



Q= / dV(V-v) 3 f(x,V;t) = 0. 



(9) 



In particular, the last expression Q implies that the 
waterbag distribution yields an exact closure of the hy- 
drodynamic equations in the form 



d t v + vd x v + v%d x (g 2 ) = d x (p, 



where 



f 



(11) 



(12) 



is the normalized thermal speed-squared. The macro- 
scopic fluid equations (fTOjl and (TTTj) must be supple- 
mented with the equation for the amplitude of the field 
mode ([3]), written in the form 

a-^+ila = i(ge- ix ), (13) 

where now the operator averaging specified by Eq. (|4]) 
involves integration on the normalized spatial variable x 
only, and the potential tp is given by expression as 
before. 

We further scale the hydrodynamic and field variables 
according to the expressions 

Q = 0o + ( R, v — v o + tV, a = ea , ^ = e^i (14) 

where e is a formal small parameter, and will be set 
equal to one at the end of the calculations. Furthermore, 
go = const and vq = const represent the stationary so- 
lution of Eqs. (|T0| and (jlll) . provided the stationary 
field amplitude ao = 0. The basic macroscopic fluid and 
electrostatic field equations can be rewritten as 

d t R + Q d x V + v d x R = -ed x (RV), (15) 



/ V 2 

8 t V + v d x V + 2 Qa v 2 d x R - d x <S> = -ed x ( — + v 2 T R 2 



a—+ila = i(Re tx ), 
dt 



(16) 
(17) 



$ = a(t)e K +a*{t)e- lx . (18) 

To simplify the above system of equations, we perform a 
Galilean transformation specified by 

z = x- v t, = a{ty v °\ (19) 

and cast our basic system of equations in the form 

d t R+g d z V = -ed z (RV), (20) 

d t V + 2g v 2 d z R - = -ed, + v 2 rA , (21) 



a-^+iCV = i(Re- lz ), 
dt x ' 



(22) 



$ = ¥(t)e 4 * + **(t) e - lz , C = l-av . (23) 

Eliminating V from the left-hand-sides of Eqs. (|2^|) and 
([2~T|) . we arrive at the basic system 



d;R-\ 2 d*R-Q <f> = -edtd z {RV)+ee d* 



V 2 



vlR 2 



d t g + d x (gv) = 0, 



(10) 



(24) 

o-St* + i&f! = i(Re~ iz ), X 2 = 1q\v%, (25) 

for the subsequent analysis using the RG approach. 
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III. RENORMALIZATION GROUP REDUCTION OF 
THE MACROSCOPIC EQUATIONS 

Following the standard procedure of the RG method^, 
we represent Q(z;t) as a perturbation expansion 



(26) 



n=0 



where Q = (R, V, $) represents all hydrodynamic and 
field variables. To zero order, the perturbation equations 
and (1231) can be written as 



d 2 R - \ 2 d 2 z R = Qo{*oe lz + %e- lz ), (27) 



ad t ^ + iC^ Q = i(R e~ lz ). 



(28) 



Since the equation for Rq contains only the first har- 
monic with regard to the spatial variable z [proportional 
to exp(±iz)], its solution is sought in the form 

Ro(z;t) = T(t)e" + J"*(t)e~ lz , (29) 

where the function J 7 satisfies the equation 

(d 2 + A 2 ) (ad t + iC)F - iqqF = 0. (30) 

Taking the latter into account, we can write 



Ro(z;t) = ^2 A„ 



ioj m t-\-iz 



(31) 



where the sum spans over all roots of the characteristic 
equation 



(uj 2 - X 2 )(auj + C) + g = 0. 



(32) 



In general, the characteristic equation has three roots. 
However, as it will become clear from the subsequent ex- 
position, in a number of cases depending on the particu- 
lar values of physical parameters one of the roots is real, 
while the other two are complex conjugate. Note also 
that the arbitrary (to this end) amplitudes A m are con- 
stants with respect to the spatial z and time t variables. 
Using Eqs. (JUJ and (2EJ), we find 



Vo = — - y^ y uj m A n .< 



go 



go 



*0 = -J2( X2 -^) A ^ mt - 

^0 _ 



(33) 
(34) 



It is worthwhile mentioning that in general, the solution 
of Eq. (|2"T|) would contain all other harmonics starting 
with the second one [proportional to exp(±mz), where 
n is integer]. They are "sound" waves, propagating at 
speed A and are decoupled from the electric field am- 
plitude "J. Their amplitudes are vanishing in the zero 



order approximation. This assertion will be verified nu- 
merically and will be discussed in more detail in Section 
V. 

In first order the basic equations (j2"4")l and (|2"5|) can be 

expressed as 

d 2 t Ri - \ 2 d 2 Rx - go®! = 



IVfi A A p*(Wm+«„)t+2i 



+ c.c, 



adt^i+iC^i =i{Rie- iz ). 



(35) 



(36) 



Here and afterwards, the symbol "c.c" stands for the 
complementary complex conjugate counterpart, while 
G mn is a symmetric matrix given by the expression 

G mn = {oJ m + w„) 2 + 2(oj m uj n + A 2 ) . (37) 

It can be verified in a straightforward manner that the 
general solution of the first order equations is 

Ri= B< + V"<* +Xt > + Yl B^e^ z -^ 

!/#±l V^±l 



IVF A A p*("m+U„)t+2 



go 



where 



(0J m + LO n f + 2(u! m UJ n + A 2 ) 



-4A 2 



(38) 



(39) 



and Bu are the arbitrary constant amplitudes of the 
solution to the homogeneous part of Eq. (|3"5j) . The sum- 
mation over v runs from — oo to oo excluding v = ±1. 
The first harmonic [proportional to exp(±iz)] describes 
the coupling between the fluid dynamic and the electric 
fields and has been taken into account in zero order. For 
the first order current velocity V\ we obtain 



Vi 



go 



J2 B(+> e ^ z+At > - £ Bi~ 



2 /J Vmji^mAi^ 

^~0 Tyi ri 



{0Jrn J rbJ n )t-\-'2iz 



+ c.c, (40) 



where 



u m w n + A 2 + 2X 2 F mn 



(41) 



Note that in first order the electric field vanishes — 0. 

The final step in our perturbative procedure is to ob- 
tain the secular second-order solution. Retaining terms 
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giving rise to secular contributions in the second-order 
solution, we can express the constitutive equations as 
follows 

d 2 R 2 - X 2 d 2 z R 2 - g <S> 2 = 



i V A nlm A* m A n A l e l ^ + ^~^ t+tz + c.c, (42) 

m , n , /. 



ad t * 2 +iCy 2 = i(R 2 e- iz ), 



(43) 



where the coupling coefficients A mn i is given by the ex- 
pression 



A-nlm = LO n LOl + A + 3A F n l 



+ (w n -u* m )u)* m F nh 



(44) 



Taking into account terms driving self- and cross- 
modulation, it is straightforward to verify that the secu- 
lar second-order solution can be expressed as 



R.2 = 

where 

T 
and 



— T A 



Qq 



e 4- c.c, 



(45) 



a(u> m + 2ij n ) + C 



frmi — A mnn , (46) 



7 m = Im(w m ), An = A m e lmt . (47) 

Moreover, the quantity V gm is the group velocity defined 
as follows. Let us introduce the dispersion function^ 3 - 

V{k,u) = (\ 2 k 2 - uj 2 ){auj + C) - g , (48) 

corresponding to a general solution of Eq. (|24[) propor- 
tional to exp (iiot + ikz), where k is the wave number. 
Obviously, the dispersion equation for k = 1 that is, 
2?(l,u;) = reduces to the characteristic equation (1321) . 
The group velocity is given by the expression^ 



V m 



dk 



dV ( dV 



k=l 

or in an explicit form 

^077 



dk \duj n 



2g a X 2 



k=l 



(49) 



(50) 



The last step is to collect all terms corresponding to in- 
creasing orders, which contribute to R(z; t) and perform 
a resummation such as to absorb secular terms (propor- 
tional to various powers of the time variable t) present 
in i?2- Since this approach is standard^ 3 -, we omit details 
here. 



IV. THE RENORMALIZATION GROUP EQUATION 

Following the procedure^ 3 - of the RG method, we fi- 
nally obtain the desired RG equation 



idtAm 



V, 



gin 



2\2 ^rnnAr, 



Ay 



(51) 



where now A m denotes the renormalized complex ampli- 
tude of the type (I47|) . In terms of the renormalized wave 
envelopes up to first order in the formal expansion pa- 
rameter, the macroscopic density g(z; t) can be expressed 
as 

Q(z; t) = g + J2 An (t)e mc ^ t+lz 



Re(cj m )t— iz 



(52) 



Similar expressions hold for the current velocity V [com- 
pare with Eq. (|33|) ] and for the electrostatic potential ^ 
[compare with Eq. (|34p]. 

It was mentioned in the previous section that since the 
characteristic equation (j3"!2)) is an algebraic equation of 
third order, it possesses three roots in general. Thus, 
we can distinguish the following two cases of physical 
interest. 

All three roots lo\, L02 and uj 3 are real. As it will be 
shown in the next Section, at certain value of the param- 
eter Vq they satisfy the resonance condition 



u)\ + lo 2 — 2oj 3 + 5. 



(53) 



Apart from the self and cross modulation terms in Eq. 
(|5ip. additional terms describing wave mixing will be 
present. Therefore, in the most general case the renor- 
malization group equation can be written as 



id t Ai 



V91 
g 2 \i 



(thIA^ 2 + r 12 \A2\ 2 + t 13 \a 3 \ 2 )a 1 



^ [u 2 + A 2 + (3A 2 + u^F^AlAle-^, 



On 



(54) 



id t A 2 = -^(t 21 \Ai\ 2 + ?22\A 2 \ 2 + T 23 \Az\ 2 )A2 
g A \ ) 



V, 



a 2 r, ,2 



8o' 



[uj 2 + A 2 + (3 A 2 + u) 1 W2)F^]AlA : 



(55) 



id t A 3 = -^2 ( r 3i I A I' + r 32 1 A | 2 + T 33 | A 1 2 ) A 3 



g 2 o\ 



c>r< 



2^( Wl w 2 + A 2 + SA 2 ^ + u; 2 f\ 3 ) AAA 

nA 



(56) 
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The terms proportional to T nn describe the effects of 
phase self-modulation, whereas those proportional to 
Tmn (with m ^ n) are responsible for phase cross- 
modulation. Both result in an intensity-dependent, non- 
linear frequency shift of the eigenfrequencies w m , an effect 
similar to the nonlinear wave-number shift in nonlinear 
optics^. In contrast, the last terms in each of the equa- 
tions above describe four-wave mixing effects, responsible 
for energy exchange between the three modes. These ef- 
fects play a role only when the phase mismatch parameter 
S is close to zero. 

One real uj\ root and two w 2 and w 2 complex conjugate 
roots. In this case the renormalization group equation 
flSTl) takes the form 



■a a ^ffl 

l9tM = ^ 



id t A 2 = - 



</2 



r 11 |^ 1 | 2 + r 12 (i + 2z 7l2 ) 



r 2 i|A| +r 22 (i + 2i 722 ) 



A±, 
(57) 

(58) 

where the quantities 712 and 7 22 are given by the expres- 
sions 



712 = 



cr 72 



<7Wl + C 



722 = 



072 



au)2 + C 



(59) 



An important comment is now in order. The equations 
for the slowly varying wave amplitudes describing the 
coupling between the hydrodynamic and electric fields 
were obtained in third order in the formal small param- 
eter. Note that up to this order the "sound" waves with 
amplitudes Bu^ remain uncoupled with the basic modes 
with amplitudes A m ■ Such coupling (similar to the wave 
mixing of the three eigenmodes) may take place provided 
a resonance condition of the form w m + oj n = ±2A is ful- 
filled, where 10 m and 10 n are any two roots of the char- 
acteristic equation ([52")) . Note also that the amplitudes 
of the " sound" waves remain undefined to third order in 
the formal expansion parameter. To obtain the corre- 
sponding amplitude equations for these modes, one has 
to go beyond third order, which is out of the scope of the 
present paper. 



V. NUMERICAL SIMULATIONS 

In this Section we present a comparison between the 
numerical solution of the system of hydrodynamic equa- 
tions (fTU)) , (jlip supplemented with the equation for the 
single wave amplitude (fT3|) . and the RG solution obtained 
in the previous section. Since the nondimensionalization 
of the initial phase space variables allows an arbitrary 
velocity scale, it is convenient to define the latter such 
that 



-14, 
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)l=0.425/- 


7.038 


Im(Q> 2 


=-Im(d> 3 ) 


=0.4248 





















2 4 6 

t 



10 12 



FIG. 1. Continuous line: time-dependence of the amplitude 
a(t) at I = — 1 and a — — 1 (one real and two complex roots of 
the characteristic equation) obtained by numerical integration 
of Eqs. ([10]), JTTJ) and ([T3)) . The asymptotic behavior of a(t) 
is given by the dashed line. Its slope is the absolute value of 
the imaginary part of the complex roots. 



In order to assess the accuracy of the model equations 
([54)) - (|56|) the predictions of the latter are benchmarked 
to the numerical solution of the original system of Eqs. 
(fT0|) . (1TT1) and (fT3|) . To achieve this an initial condition 
of the form g(x; t = 0) = 1 + e cos(x) with e in the range 
0.001 . . .0.01 and v(x;t = 0) = v n = const(x) has been 
used. Equations ([3"Tj) . (|33l) and ([54"f relate the initial 
distributions of g(x;t = 0), v(x;t = 0) and fy(t = 0) to 
the three slowly varying amplitudes A\(t = 0), *4 2 (t = 0) 



and Az{t — 0). Varying e in the range 10 



-3 



10" 



(60) 



cor- 
responds to a transition from a linear to a weakly nonlin- 
ear system dynamics. Periodic boundary conditions have 
been implemented. 

The results shown in Figures[T]and[2]pertain to the lin- 
ear dynamics regime. It can be shown that in the case, 
where / = —1 and a = — 1 (Fig. 1) the characteristic 
equation (f3"2"j) has one real and two complex conjugate 
roots. From Eq. (|34|) it is clear that the mode corre- 
sponding to the root with negative imaginary part will 
become dominant for sufficiently late times t. This ob- 
servation is in full agreement with Fig. 2. In addition, 
the slope of the asymptotic line matches the imaginary 
part of the complex roots with an accuracy of the order 
of 10~ 4 . In contrast, when I = 1 and a = 1 the character- 
istic equation possesses three real roots and these corre- 
spond to the three eigenfrequencies of the system lo\ , w 2 
and W3. By Fourier-analyzing the numerical data for the 
amplitude of the electrostatic potential a(t), the values 
of these frequencies can be obtained and compared to the 
roots of the characteristic equation. As Fig. [2] shows the 
two results are in excellent agreement. 

Increasing e from 10~ 3 to 10~ 2 effectively switches 
on the nonlinearity. It can be shown that the phase- 
matching (resonance) condition (|53l) is satisfied at v = 
Vpm ~ 0.49. In order to assess the efficiency of the four- 
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o=l, 1=1, v =0, a(0)=0, 8=0.001 
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FIG. 2. Spectrum a(cj) of the amplitude a(t) in the case of 
1 = 1 and cr = 1 obtained from the numerical solution of Eqs. 
(|10|l . and (j 13 P - The three peaks correspond to the three 
real roots of the characteristic equation. Note the absence of 
"sound" waves at frequencies uj = ±\/3. 
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FIG. 3. Fourier transform v(uu,k = 2) (k is the harmonic 
number) of the current velocity v(z;t). The dashed line is the 
spectrum obtained by a fast Fourier transform. In addition 
"slow" Fourier transform is performed in the vicinity of the 
peaks in order to improve the frequency resolution and the 
continuous line is the result. The "*" -symbols represent the 
analytical result, given by Eq. (|40p . The values of the eigen- 
frequencies are U\ = 1.589, u>2 = —0.173 and U3 = —1.925. 
The two components at frequencies — 2A and 2A are "sound" 
waves (harmonic number k = 2), propagating at speed A. 



wave mixing (FWM) process, Eqs. ([54)) - (|56|) have been 
solved for several values of vq in the vicinity of v pm . It 
was found that significant energy exchange between the 
modes occurs on time scales of the order of 10 4 . Solv- 
ing the original system of hydrodynamic equations ([10j) , 
(ITT|) and IT51) on such a timescale is not practical and 
given the small efficiency of the FWM process it can be 
neglected. 

Figures |3] and U compare the numerical data for the 
current velocity v(z;t) with the analytical result, given 



o=l, 1=1, v =0.49, a(0)=0.0055, 8=0.01 
x10 



-0.36 -0.35 -0.34 -0.33 -0.32 

CO 

FIG. 4. The same as in Fig. [3] but the frequency region near 
2oi2 and ui\ + u>3 is shown in more detail. 



by Eq. (|40)) . Firstly, the quantity v(u,k = 2) has 
been computed by fast-Fourier transforming the quantity 
(v(z; t) cxp(— 2iz)) and this results in locating the peaks 
of the spectrum (dashed line in Figs. |3]and|2]). In order to 
increase the accuracy " slow" Fourier transform has been 
performed in the vicinity of the peaks (continuous lines) . 
Finally, the "*" -symbols represent the amplitudes of the 
components at frequencies u m + u> n computed from Eq. 
(|40|) . As can be seen, the analytical and numerical re- 
sults are in very good agreement. Note, that besides the 
six frequency components at frequencies u> m + ui n there 
are two additional peaks at frequencies u> = ±2A. They 
represent "sound" waves, solutions of the Eqs. (|T0|) . (fTTj) 
and (|13|) with tp — and harmonic number k = 2 prop- 
agating at the sound velocity A. Since these waves are 
not in resonance with any of the six components at fre- 
quencies U) m + u> n (±2A u) m + U) n ) and their amplitudes 
remain much smaller than A\ , A2 and ^3 , their effect on 
the evolution of the latter can be neglected. 

Figure O compares the analytical and numerical results 
for the density g(x;t). The curve labeled "Linear" is 
obtained from Eq. (|31[) with A m — const. The curve 
with label " Nonlinear" in addition takes into account the 
" slow" variation of A m , resulting from the nonlinear self- 
and cross-modulation effects described by Eqs. (|54[) - 
(|56|) . As can be seen taking into account the nonlinear 
effect improves the accuracy of the analytical solution. 
This has been verified by calculating the square difference 
between the numerical and the "linear" and "nonlinear" 
analytical results, respectively. 



VI. CONCLUDING REMARKS 

Using a class of waterbag phase space density distribu- 
tions, which is an exact solutions to the Vlasov equation, 
we have cast the single-wave model equations to a fluid 
dynamic form with nonzero thermal velocity. Since the 
continuity and momentum balance equations comprise 
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FIG. 5. "Numerical": density distribution g(x;t = 600) ob- 
tained from the numerical solution of Eqs. (|10l) . (JTTJ and (I13|) ; 
"Linear": Eq. (I31[) with A m = const. The curve marked 
"Nonlinear" is obtained by taking into account the nonlinear 
frequency shift of Am, Eqs. (JS3J - ((56)) . 



an exact hydrodynamic closure, the hydrodynamic rep- 
resentation is fully equivalent to the original system com- 
prising the Vlasov equation and the evolution equation 
for the single electrostatic field mode. 

Based on the renormalization group method, a sys- 
tem of coupled nonlinear equations for the slowly varying 
amplitudes of interacting plasma density waves has been 
derived. Depending on the solution of the characteristic 
equation, the system mentioned above takes the form of 
either three coupled nonlinear equations describing the 
process of three-wave interaction in the case, where all 
three roots are real, or two coupled equations with com- 
plex coefficients if the characteristic equation possesses 
one real and two complex conjugate roots. 

In linear dynamic mode (and three real roots of the 
characteristic equation) , the values of the three eigenfre- 
quencies u>x, 0J2 and W3 computed numerically have been 
found in excellent agreement with the solutions of the 
characteristic equation (j32j). Similarly, in the case of a 
single real root, the asymptotic behavior of the electro- 
static potential is in full accordance with the value of 
imaginary part of the complex roots. In the weakly non- 
linear regime, the spectrum of the nonlinear correction 
to the current velocity closely matches the analytical so- 
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lution. Finally, direct comparison between the analytical 
expressions for the density g(x; t), with and without non- 
linear effects, and its numerical solution via numerical 
integration of the hydrodynamic equations shows that 
taking the nonlinear effects into account improves the 
agreement between the analytical and numerical results. 

It is worthwhile to mention that similar equations have 
been obtained earlier by Janssen and Rasmusseni. Al- 
though in a different physical context, amplitude equa- 
tions similar to the ones derived in the previous section 
have been reported by B. Bruhn et a l 15 ' 16 . Recently, 
numerical simulations of a discrete self-consistent wave- 
particle model describing an ensemble of particles cou- 
pled to a monochromatic wave have revealed a new phe- 
nomenon, a phase transition associated with the Landau 
damping regime. This latter finding adds particular proof 
to the result presented here. 

An interesting continuation of the present study would 
be the consideration of large-scale hydrodynamic fluctu- 
ations and their influence on the dynamics of perturba- 
tions. This we plan to complete in a future publication. 
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